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, The analytical probability distribution of the quasi-2D (and purely 2D) ideal and interacting Bose 

' gas are investigated by using a canonical ensemble approach. Using the analytical probability distri- 

I bution of the condensate, the statistical properties such as the mean occupation number and particle 

> I ■ number fluctuations of the condensate are calculated. Researches show that there is a continuous 

' crossover of the statistical properties from a quasi-2D to a purely 2D ideal or interacting gases. 

Different from the case of a 3D Bose gas, the interaction between atoms changes in a deep way the 
nature of the particle number fluctuations. 
00 ' hongweixiong@hotmail.com 
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I. INTRODUCTION 



The remarkable observations of Bose-Einstein condensation (BEG) in atomic vapors have opened up new 

avenues of research into the physical properties and nature of Bose-condensed systems [pppT Recently BECs have 
been realized in quasi-one and quasi- two dimensions [^,^, where new phenomena such as quasicondensates with a 
fluctuating phase |p|-|ri| and a Tonks gas of impenetrable bosons ||l^,|l2|,|l3| may be observed. After the realization of 
^ , BEG, the statistical properties of the system (especially the particle number fluctuations of the condensate) have been 
Q ' studied in deep for three-dimensional (3D) ideal Bose gases confined in a box |l4|jr^ ], and in the presence of a harmonic 
O ] trap |p^-p2| . Recently the question of how interatomic interactions affect the particle number fluctuations has also 
been the object of several theoretical investigations [p3|-p0[. Recall that BEGs can not occur in one-dimensional (ID) 
' and two-dimensional (2D) uniform Bose gases at finite temperature, because in that case thermal fluctuations would 
^ . destabilize the condensate. The realization of the trapped BEGs in various dimensions makes the behavior of particle 

■ number fluctuations a very interesting problem, especially for ID and 2D Bose gases in a harmonic trap. 

, It is well known that the theory based on a grand-canonical ensemble gives unphysically large particle number 
' fluctuations of condensate, in particular when the temperature approaches zero. It is therefore necessary to use a 
\ canonical (or microcanonical) ensemble to investigate the statistical properties of the trapped Bose gas. Recently a 

■ simple method has been developed p9|-|3l|] to calculate the probability distribution of the condensate within canonical 
] ensembles. As soon as the probability distribution of the condensate is known, the statistical properties of the 

"Y"; ■ system can be obtained directly. With this method the mean ground state occupation number and particle number 
d ' fluctuations for the 3D have been considered for interacting gas trapped in a box |^ and a harmonic potential 
^ . The role of dimensionality on the statistical properties of the ideal gas was also discussed recently through 

I the calculations of the probability distribution for different spatial dimensions and various confining potentials. The 
' ^ , purpose of the present work is an attempt to investigate the probability distribution of quasi-2D (and purely 2D) 
ideal and interacting gases confined in a harmonic trap. Using the probability distribution function obtained, we can 
calculate the statistical properties of the system such as the mean occupation number and particle number fiuctuations 
of condensate in a simple way. 

The paper is organized as follows. In Sec. II we give a simple introduction of our approach to calculate the statistical 
properties of the Bose gas within the canonical ensemble. In Sees. Ill and IV we investigate the statistical properties 
^ of quasi-2D (and purely 2D) ideal and interacting Bose gases, respectively. The last section (Sec. V) contains a 
5^ , discussion and summary of our results. 

II. PROBABILITY DISTRIBUTION WITHIN CANONICAL ENSEMBLE 

The statistical method used in this work follows from that used for the 3D dilute interacting gases developed in 
Ref. [ p9|j30| ] , where a canonical ensemble is used to discuss the statistical properties of the system. We start from the 
canonical partition function of N trapped interacting bosons, given by 

Z{N)^ exp[-/3(i;o + I]„^oA^„£„ + i;,„0], (1) 

E„JV„=Ar 



1 



where /3 = 1/kBT, Nn and e„ are occupation number and energy level of the state n, respectively. Eq is the energy 
of the condensate of the system. In the case of non-interacting Bose gas, Eq = NqEq; for the case of interacting Bose 
gas, however, Eq should be regarded as the interaction energy of the condensate. In the above equation, Eint is the 
interaction energy between the condensate and the normal gas, which takes the form 



E^nt ^2g J no (r) ut (r) d^r + g J n\ (r) d^r. 



(2) 



where no (r) and ut (r) are the density distribution of the condensate and normal gas, respectively. In the above 
equation, g represents the coupling constant. 

The probability distribution of a condensate is defined as the probability to find A^o atoms in the condensate. 
Through a developed saddle-point method |p9| , ^ one can obtain the probability distribution of the condensate, 
which is found to be 



G„(iV,Aro) = A„exp 



a(iV,iVo)diVo 



(3) 



where An and iVp are normalization constant and the most probability value of the condensate, respectively. The 
quantity a (TV, TVq) can be obtained based on the following equation: 



exp exp[-/3g^ {Eq + E.nt) - a {N, No)] - 1 



By setting a {N,No) = in Eq. (jj) one gets the most probability value Nq, 

1 



exp [/3e„] exp[-/3g^ {Eq + Ei^t)] - 1 ' 



(4) 



(5) 



Note that Nq is exactly the mean ground state occupation number within a grand canonical ensemble, when the 
lowest order perturbation theory is used to discuss the role of atom-atom interactions (see Ref. [^,^,^ ) . Although 
a {N, No) can be obtained directly from Eq. (^, it is more convenient to calculate it through the difference between 
No and Nq, i.e., 

No-N'^^J2 



exp [Ps-n] exp[-/3g^ {Eo + E.m)] - 1 



J„ exp [/?£„] exp[-pj^ {Eo + Ei^t) - a {N, No)] - 1 

Once we know £„, Eo and Eint of the system, it is straightforward to obtain a{N,No) from Eq. (^. Then by Eq. 
(^ one can obtain the probability distribution of the condensate directly. 

As soon as we know Gn {N, No), the statistical properties of the system can be clearly described. From the 
probability distribution function given by Eq. ^ one obtains (A^o) and (^6'^No) within the canonical ensemble: 

ENo=oNoGn{N,No) 
^"""^ ^ El.oGAN,No) ' 



ENo=oNiGn {N,No) 

ENo=oGn{N,No) 



(S'No) = 



ENo=oNoGn {N,No) 

ENo=oGniN,No) J 



(8) 



We see that in this approach, the calculation of a {N, No) by using Eq. (^ plays a crucial role to discuss the statistical 
properties of the condensate. 
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III. IDEAL BOSE GASES IN A QUASI-2D HARMONIC TRAP 



Now we apply the formulation introduced in the last section to calculate the thermodynamical quantities such as 
the particle number fluctuations of the condensate for an ideal 2D Bose gas. We first consider the statistical properties 
of a purely 2D Bose gas confined in axially symmetric traps. The purely 2D Bose gas can be regarded as the limit of 
a trapped 3D Bose gas with the trapping frequency lo^ in the z direction approaching infinity. 

The single-particle energy level of the purely 2D Bose gas in an axially symmetric harmonic trap takes the form 



En = {nx +ny + l) hujj_, 

where uj± — ui^ — ojy is the trapping frequency in the radial direction. From Eqs. and (|^), one gets 

1 1 



exp [l3 {r 



j) hLu±] — 1 exp [13 {hx + ny) %tjj^_] exp[— a (TV, Nq)] — 1 



(9) 



(10) 



Assuming \a {N, No)\ << 1, we obtain the result of a {N, Nq) 

a {N, No) ^ 

where 



No - < 



^2D — 



E 



exp [P {rix + Uy) fiLin} 
{exp [(3 (rix + Hy) hLLi±] - 1}^ 



(11) 



(12) 



It is easy to find that Eq. (|ll| ) coincides with the assumption |a(iV, A^o)| << 1- Using Eq. (H), the probability 
distribution of the purely 2D ideal gas is then given by 



G2D = Md exp 



{No-N^oY 



2S 



2D 



(13) 



where A2D is a normalization constant. 

In addition, from Eqs. (|^) and (|9|), we obtain the most probable value Nq as 



N- 



-Nt^ 



tN^'^ In TV, 



(14) 



where we have introduced a reduced temperature t ~ T/T^, with ~ {N/C,{2))^^'^ TiLj^/kB being the critical 
temperature of the system in the thermodynamic limit. 

Using the formulas (|), (|), (|l3|) and (0), we can obtain the mean ground state occupation number and particle 
number fiuctuations of the condensate. In particular, below the critical temperature (i.e., in the case of A^q >> 1), 
the analytical result of the particle number fiuctuations is given by 



{S^No) 



2D 



-^2D 



Ne 
W) 



■In 



Ne 
W) 



(15) 



When obtaining Eq. (p^) we have used the formula (^S^Nq) =a/2 for the particle number fluctuations when the 



«-iVo)Va 



and N'q » 1. 



A 



probability distribution is a Gaussian function with the form G„ = yl„ exp 

simple derivation for this result has been given in the Appendix of Ref . |]3l| . 

We see that the particle number fiuctuations of the harmonically trapped 2D gas exhibit a very weakly anomalous 
behavior, as it is controlled by the factor In {^Nt^ /C^ (2)) . Shown in Fig. 1 is the numerical result of 5Nq = yJTPNo)^ 
(solid line) for N = 10"^ purely 2D ideal bosons confined in an axially symmetric trap. The numerical result of 
(A'o) /N (solid line) for N = 10"^ purely 2D ideal bosons is demonstrated in Fig. 2. As a comparison, (A^o) /N within 
a grand-canonical ensemble (or Nq within the canonical ensemble) is also shown in the figure. It is easy to show that 
Nq coincides with (A^o) in the case of iV^ >> 1. However, there is an obvious difference between Nq and (iVo) near 
the critical temperature when the total number of particles iV is smaller than 5 x 10"^. 
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exp [j3 {rix + Uy) Tiuja] exp[— a {N, Nq)] — 1 



oc 

E 

1 Tlx j fly —0 



E 



For convenience we separate off the state n 
By setting a{N,No) = in Eq. ([F 



exp [P [n^ + rtj,) Tmj±\ exp [iSrizhuJz] exp[— a (TV, Nq)] — 1 

from the state Uz ^ in the sum on the right hand side of Eq. (|1 



(17) 



we get Nq-. 



N - 



exp [f3 {fix + Tiy) huj^] — 1 



CX) oo 

E E 



exp [P {rix + Hj/) ?ia;j_] exp [/Suzhajz 
From Eqs. (|l|) and it is straightforward to get the finial result of a {N, Nq) 

a (N, No) 



No 



^2D 



e 



2D 



where is the constant determined by Eq. (|12|), and 

^ exp[f3 {rix + ny)huj±]exp[f3nzUujz 



e 



„,=i„,.„^=o{6xp [/3 ("x + riy) hLj±] exp [/Jn^Titj^] - 1}^ 
For the quasi-2D gas, we have ksT < fiujz. In this situation, Q2D can be approximated as: 



e 



2D 



C(2) 



E 

.ji^=i 

The probabihty distribution of the quasi-2D gas is then 

Gq-2D ^ Aq^2D exp 



exp {-PuzfiUz) 



(No ~ NEf 

2{E2D+e2D) 



(18) 



(19) 



(20) 



(21) 



(22) 



where Aq^2D is a normahzation constant. In addition, from Eq. (|l 
as: 

N-Nt'^ 



the most probabihty value can be approximated 



N'o 



'tN^/^\nN -e2D- (23) 
represents a modification between a quasi-2D and a purely 2D 



The last term 82D on the right hand side of Eq. (| 

Bose gas. 

From the formulas (|), (|), (H) and (||), we can obtain the mean ground state occupation number and particle 
number fluctuations of the condensate for the quasi-2D gas. Below the critical temperature, the particle number 
fluctuations of the condensate for the quasi-2D gas is given by 

{S'No)^_,^ 



^2D 



e 



215- 



(24) 

represents a modification 



Similar to the most probability value, the last term Q2D on the right hand side of Eq. (^ 

between the quasi-2D and the purely 2D Bose gas. 

In Fig. 1, the numerical results of 6 No for the quasi-2D ideal gas are shown for z = 1 (dashed line) and z = 2 (dotted 
line), respectively. The continuous crossover of 6N0 from the quasi-2D to the purely 2D gases is clearly illustrated in 
the figure. Below the critical temperature, we see that 6N0 of the quasi-2D gas is slightly larger than the result of 
the purely 2D gas because confinement has the effect of reducing the particle number fluctuations. As a comparison, 
the dot-dashed line in Fig. 1 displays SNo oi N = 10^ ideal bosons confined in a 3D istropic harmonic trap, where 
the oscillator frequency is assumed to be equal to the transverse frequency lu± of 2D Bose gas. The arrow marks 
the critical temperature of the 3D Bose gas in the thermodynamic limit. Shown in Fig. 2 are the numerical results 
of (iVo) /N for the quasi-2D ideal gas with z = 1 (dashed line) and z = 2 (dotted line). There is also a continuous 
crossover of (iVo) /N from quasi-2D to purely 2D gases. In addition, (No) /N for the quasi-2D gas is slightly smaller 
than the result of the purely 2D gas. 
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IV. INTERACTING BOSE GASES IN A QUASI-2D HARMONIC TRAP 



In this section we turn to discuss the statistical properties of a quasi-2D interacting Bose gas in an axiaUy symmetric 
trap. The interaction between atoms can be described by a single coupling constant g « 2\/27r?i^as/ {ml^) which 
is fixed by a s -wave scattering length and the oscillator length 1^ = (Ti/mWz)^/^ in the z direction. Note that the 
coupling constant in the quasi-2D interacting gas is related to the trapping frequency oJz in the z direction because it 
is the trapping frequency lo^ that leads to the crossover from a 3D to a quasi-2D gas. This unique property makes the 
statistical properties of the quasi-2D interacting gas a very interesting problem. By the investigation of the statistical 
properties of the quasi-2D interacting gas, one can get a clear description of the role of the trapping frequency For 
a weakly repulsive interacting Bose-condensed gas, the correlation length h/^mnog (uq is the density distribution of 
the condensate) should greatly exceed the mean interparticle separation 1/ ^27rno- In this situation, one can introduce 
a small parameter 6 = mg/2'K?)^ « 1 for the weakly interacting gas 

In Thomas- Fermi regime the chemical potential of the quasi-2D gas is given by fi{No,T) — {Nomg/ny/^Ljj^. 
In addition, there is an important relation between Eq and fi{No,T): /i(iVo,T) — dEo/ ONq. It is convenient to 
introduce a dimcnsionless parameter 77 (A^o, T) defined by 



It can be rewritten as: 



r;(7Vo,r)=r]o (!-<') A, 



(25) 



(26) 



where t is the reduced temperature introduced before and 770 = a/2(' (2) 6 = \.&2{as/lzY^'^ is a small parameter for 
the weakly interacting gas. Note that 770 is irrelative to the number of the condensed atoms. For the low-dimensional 
experiment |Q on Bose-condensed systems a simple calculation gives 770 = 0.1. We will see in the following that the 
role of atom-atom interactions on the thermodynamic properties of the system is determined by the parameter 779. 

For the quasi-2D interacting gas, with the lowest order approximation, the occupation number N,^ of the exited 
state n 7^ is given by 



J huji] exp[Pnzhujz-r] {Nq, T)] - 1 ' 



expl 



(27) 



The term rj (7Vo,T) in the above equation represents the interaction correction due to the interatomic interacion in 
the condensate. Different from the ideal gas, we see from Eq. ( |27| ) that the occurrence of the quasi-2D interacting 
gas should satisfy the condition tkoz >n{No,T) » huj± (see also Ref. Q). Note that the occurrence of a quasi-2D 
ideal gas should satisfy the condition hujz >kBT >> fiuj±. 

Omitting the interaction between condensed and normal gases, from Eq. (^, we get 



No = N- J2 



1 



exp [P {n^ + uy) Huja] exp [-77 (iVo, T)] exp[-a (iV, A^o) 



00 00 

E E 



n^—l Tij. .n^— 



exp [P {ux + Uy) hu!±] exp [jSrizhuJz - i] (No, T)] exp[-a {N, No)] - 1 
By setting a{N,No)^0 in Eq. (|2|), is given by 

1 



(28) 



exp [f3 [iij, + Uy) Tiuj±\ exp [-rj (A^o, r)] - 1 



oo oo 



n,— 1 rij. .n-u—0 



E E 

From Eqs. (M) and (M), we get 



exp [P {ux + Uy) fiuji} exp [PuzfiUz - rj {Nq, T)] - 1 



(29) 
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(30) 



where 



T2D — 



exp [/3 {n^ + Uy) huj± - 77 {Nq, T)] - 1 



and 



exp [/? ( Tlx 



^ mod ^ ^ ^ 



c« oo 



n, — 1 n^., 71^,-0 



exp [/3 {n.j; + Uy) hu;±] exp [l3n^huj^ - rj {Nq, T)] - 1 



(31) 



exp [/? {ux + riy) hLu±] exp [f3nzhuJz - V {No, T)] exp[-Q; {N, Nq)] - 1 



(32) 



A. Pure two-dimensional interacting Bose gases 

Firstly we discuss a purely 2D interacting gas where Tmod can be omitted because f3hujz 
Thus for the purely 2D interacting gas, from Eq. (pOj), one gets 



in this situation. 



W) 



[t? (iVo, T) In 77 (iVo, T) - r; «, T) In 77 «, T)] + a In 7/ (TVq, T)} 



(33) 



When obtaining Eq. (p3), we have used rj {N^, T) << 1 by assuming that the interaction is weak. In addition, when 
obtaining this formula we have used the assumption a(iV, A^q) <<^(-^0'^) which may not hold in the case of an 
extremely weak interaction. However, we will see in the following that this assumption is reasonable over a wide range 
of the experimental parameters. 

Using the following Taylor expansion 



ry(7Vo,r)=r7«,T) 
one gets the result for a {N, Nq): 



171 « - No) , 



(34) 



a [N, No) 



l + mt\lnv{NS,T)\/ 2C(2)(l-i2) 



,1/2- 



Nt^\\n7j{N^,T)\/Ci2) 
Below the critical temperature, a (TV, A^o) can be approximated as: 

No-NP 



(No-NE). 



a (TV, A^o) 



Nt^\\n7^{NP,T)\/C{2)- 



(35) 



(36) 



We see that the assumption a {N, No) <<i] {N^, T) holds in the present low-dimensional Bose condensed experiment, 
where r]Q — Q.l M. The distribution function of the condensate below the critical temperature is then 



G2D = Md exp 



{N^ - N^of 



2iVi2|in^(iVo^T)|/C(2)_ 
where A2D is a normalization constant. In addition, the most probability value Nq is given by 



(37) 



7 



□ 

D D □ □ 



□ 



B. Quasi- two-dimensional interacting Bose gases 



Now we discuss the statistical properties of a quasi-2D interacting gas. For the quasi-2D interacting gas, the second 
term Tmod on the right hand side of Eq. ( ^ ) can not be omitted. In the case of huj^ >^{No,T), Tmod can be 
approximated as: 



2D, 



(40) 



where Q2D is the constant given by Eq. (^). Using the approximation of T„iod it is easy to obtain the probability 
distribution and particle number fluctuations of the quasi-2D interacting gas below the critical temperature, which is 
given by 



Gq-2D = Aq^2D GXp 



2iVi2|lnr;(iV?;,T)|/C(2) + 2e2D 



(41) 



\^ ^0)q-2D - ^(^) + ^2D. (42) 

In addition, the most probability value is given by 

= N-Nt"^ ~ Nr]o {l - f) t [l ~ In {rio (l - t^) A)] /C (2) - iA^^^^ In - 620. (43) 

Based on the formulas (0), (||), (|l|), and (||), the mean ground state occupation number and particle number 
fluctuations of the condensate are calculated numerically. In Fig. 3, the numerical results of (iVo) /N for TV = 10'^ 
quasi-2D interacting bosons with rjo = 0.1 are shown for z = 1 (dashed line) and z = 2 (dotted line), respectively. 
Shown in Fig. 4 is the numerical result of SNq with the same interaction parameter rjQ = 0.1. Similarly to 2D ideal 
gas, below the critical temperature, we see that the particle number fluctuations in a quai-2D interacting gas is slightly 
larger than the result of a purely 2D one. 



V. DISCUSSION AND CONCLUSION 



In this work, the statistical properties of pure- and quasi-2D Bose gases are investigated within the canonical 
ensemble by using the method developed recently in Refs. p9|-^. It is found that there is a continuous crossover 
of the statistical properties from quasi-2D to purely 2D gases. Different from the case of a 3D gas, the effects of 
two-body interaction and strong confinement in the z-direction change drastically the nature of the particle number 
fluctuations for a quasi-2D Bose gas. For the quasi-2D interacting gas the role of atom-atom interaction on the 
particle number fluctuations is determined by the reduced interaction parameter rjo-, which is determined not only 
by the s-wave scattering length but also by the trapping frequency uj^ in the z direction. Thus one can control the 
particle number fluctuations in the condensate by adjusting the trapping frequency in the z direction. This unique 
property makes it very promising to observe the particle number fluctuations in quasi-2D condensates. 

In the present work, the contribution to the particle number fluctuations due to collective excitations is omitted. For 
3D Bose gas, as pointed out in the role of collective excitations would be very important near zero temperature. 
By using Bogoliubov approach, the role of the collective excitations on the condensate fluctuations would be very 
interested for the case of 2D Bose gas. After the realization of low-dimensional Bose gas, the property of collective 
excitations in 2D Bose gas becomes a very interested problem. With the development of the theoretical research on 
the collective excitations of 2D Bose condensed gas, we anticipate that the role of the collective excitations on the 
condensate fluctuations would be given in the near future. 

It is obvious that the method given here can be also applied to other Bose systems, such as an elongated (quasi-lD) 
trap. The discussion on the statistical properties of the quasi-lD interacting Bose gas is beyond the scope of this 
work and will be presented elsewhere. 
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